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Abstract 

Understanding the functional role of invasive species is important for better management of ecosystems. 
Exotic willows (Salix spp.) are widespread riparian trees of temperate Australia and New Zealand. Despite be¬ 
ing classed as weeds of national significance, little is known of their effects on ecosystem function. Extensive 
growth of complex willow root mats in streams they infest has created a novel littoral habitat. Spatial and tem¬ 
poral changes of macroinvertebrates in willow root mats and bare bank habitats were examined in three central 
Victorian rivers to gain an understanding of their effects on this faunal group. Macroinvertebrate richness, 
abundance, biomass and diversity were significantly higher (p<0.05) in willow root habitats compared to open 
bank habitats. This suggests willow root mats have created a new niche for macroinvertebrate communities in 
willow-infested streams in Australia and emphasises the need for concurrent habitat enrichment programs 
when removal of exotic vegetation takes place. (7 he Victorian Naturalist 132 (4) 96-107) 
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Introduction 

Habitat diversity has been correlated with high 
macroinvertebrate diversity (Schiemer et al 
1991; Minshall 1988; Harper et al 1992; Min- 
shall and Robinson 1998) suggesting loss of 
habitat complexity due to river channelisation 
and removal of riparian vegetation, which will 
adversely affect macroinvertebrate richness, 
abundance, biomass and diversity. This suggests 
that care is required where clearing of already 
established exotic riparian vegetation is under¬ 
taken. Recognition of available habitat struc¬ 
tures or the presence of surrogate habitat struc¬ 
tures is important for appropriate management. 

There are many examples from marine and 
freshwater environments where a positive in¬ 
crease in diversity has been achieved through 
habitat manipulations. Artificial reefs and other 
habitat structures are widely used to improve 
marine habitats to improve fisheries (Russell 
1976; Talbot et al 1978; Burchmore et al 1985). 
Bell et al (1985; 1987) used artificial sea grass 
units to increase recruitment of post larval 
and juvenile fish and crustaceans to estuaries 
around Sydney, NSW. Koehn (1987) showed 
increase in fish abundance in a stretch of the 
Ovens River, in north-eastern Victoria, after it 
had been seeded with large rocks. Snags have 
been identified as important habitat for inver¬ 


tebrate and vertebrate biodiversity (Borchardt 
1993; Phillips and Kilambi 1994). Habitat en¬ 
hancement through planting of willow Salix 
sp. and Common Reed Phragmiles australis in 
off-channel bays in parts of the ITuntspill River, 
Somerset, England, resulted in a significant in¬ 
crease in abundance and diversity of fish (Lan- 
gler and Smith 2001). 

Willows were introduced to Australia during 
the early 19th century and became naturalised. 
They have spread across approximately 30 000 
km of the 68000 km river frontage in Victoria 
(Ladson et al 1997). The impacts of willows 
on stream channel morphology have been 
widely researched (c.g. Young 1980; Ladson 
et al 1997); but little research has been carried 
out to understand their effects as habitat. Wil¬ 
lows have a characteristic complex root system 
that often grows into the stream channel and 
is different from the root system of most na¬ 
tive riparian tree species. Where willows are 
dominant, their roots cover a large area of the 
littoral zone of channels. These root mats can be 
considered as a novel littoral habitat. Studies in 
the Murrumbidgee River near Wagga Wagga, 
NSW (Robertson 1993) revealed that exposed 
inundated roots formed patches of habitat pro¬ 
viding shelter, food and oviposition sites for 
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diverse aquatic macroinvertebrates while bare 
banks were inhabited by only thinly scattered 
oligochaetes and chirononiids. Boulton and 
Lloyd (1991) also found bare banks to be poor 
microhabitat in the lower Murray River com¬ 
pared to root habitats. Their conclusions were 
based on a study of only two seasons of the year. 
Jayawardana and Westbrook (2010) and Jaya 
wardana (2011) compared macro invertebrate 
communites of root habitats provided by wil¬ 
lows with those provided by native vegetation 
and found differences in species composition. 
No comprehensive study has been carried out 
to investigate the contribution of willow roots as 
a habitat for macroinvertebrate communities in 
places where other vegetation is lacking. In the 
management of invasive riparian vegetation, it 
is important to understand the habitat function 
of these introduced species and their ecosystem 
functions to aid management and habitat reme¬ 
diation. Therefore, the objective of this study 
was to investigate how macroinvertebrate com¬ 
munities use willow root mats as a habitat com¬ 
pared to open bank habitats. Macroinvertebrate 
richness, abundance, biomass, diversity and 
assemblages associated with willow root mats 
and open bank habitats were compared in three 
Central Victorian streams during summer, au¬ 
tumn, winter and spring of 2004. 

Materials and methods 
Study sites 

Birch Creek, located in the north central re¬ 
gion of Victoria, is a major tributary of the 
northward flowing Tullaroop Creek system in 
the upper catchment of the Loddon River, The 
riparian vegetation includes native Eucalyptus 
spp.. Acacia spp., Callistemon spp>, Leplasper- 
mum spp. and exotic Salix spp. (willows), with 
intermittent reed-dominated stretches. The 
dominant land uses adjacent to the creek are 
grazing and cropping. )im Crow Creek is also 
a tributary of the Loddon River, containing 
native vegetation, exotic willows and intermit¬ 
tent reed beds. It flows largely through public 
land. The Last Moorabool is a branch ol the 
Moorabool River, which originates in the Great 
Dividing Range in Central Victoria. It flows 
through native forest and agricultural areas and 
has stretches dominated by willows and native 
vegetation (Fig.l). 


Sampling 

Two sites (100 m reach) were selected for sam¬ 
pling in each of three rivers (4th to 5th order), 
each site having stretches of willows and bare 
banks. Six samples were collected from each 
habitat type. Sampling was carried out dur¬ 
ing summer, autumn, winter and spring using 
a 20 cm x 20 cm horizontal metal frame at¬ 
tached to a sampler with a 250 pm mesh net. 
The area inside the sample frame was excavated 
to 5 cm depth and the whole sample, together 
with the portion retained in the net as a result 
of site disturbance, was collected and sealed in 
a polythene bag. All samples were transported 
to the laboratory for sorting. Dissolved oxygen, 
pH, conductivity, temperature and flow also 
were recorded in each instance. In the labora¬ 
tory, samples were sieved (minimum sieve size 
250 pm) and macroinvertebrates were sepa¬ 
rated and preserved in 70% alcohol. The spe¬ 
cies present in each sample were identified to 
the lowest possible taxonomic level (Cranston 
1996; CSIRO 1999; Gooderham and Tsyrlin 
2002; Hawking and Theischinger 1999). The 
live willow roots were separated and the re¬ 
maining fraction of the sample was separated 
into coarse particulate matter fraction (CPOM) 
(>1 mm) and fine particulate matter fraction 
(FPOM) (1 mm-250 pm), air dried for 24 hrs 
and oven dried at 105°C until constant weight 
was attained. After recording the dry weight, 
each fraction was ashed at 550°C for six hours 
and ash free dry weight (AFDW) was recorded 
for each fraction. CPOM and FPOM content 
were calculated for each sample and percent¬ 
ages of CPOM and FPOM were calculated. 

Analysis 

A total of 288 samples were collected over the 
four seasons. The data were analysed using mul¬ 
tifactorial ANOVA with season, river and habi¬ 
tats within sites as main effects. Differences in 
species richness, total abundance and Shannon 
Diversity Index (SDI) between habitats were as¬ 
sessed using multifactorial ANOVA in the ‘R’ 
package (R Development Core Team 2004). To¬ 
tal abundance values were log transformed [log 
(1 +x)] before analysis to counteract the effect of 
highly dominant taxa. 

Macroinvertebrate community assemblage 
differences were assessed using Canonical 
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Fig. 1. Sampling sites of Birch Creek, Jim Creek and Moorabool River 
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Analysis of Principle Co-ordinates (CAP). 
This technique allows (i) testing for differences 
among assemblages (ii) identification of taxa 
most highly associated with the effects of habi¬ 
tats (iii) visualisation of the between and within 
group variation and (iv) a discriminant analy¬ 
sis (Anderson and Willis 2003; Anderson and 
Robinson 2003). 

CAP analysis of communities was conducted 
in two ways, first, discriminant analysis was 
conducted to discriminate between habitats 
(willow and bare banks) for each river during 
different seasons. Confusion tables ([mis] clas¬ 
sification rates) were also generated. Secondly, 
CAP was used to determine taxa associated 
with each habitat type. For that, CAP analysis 
was conducted for two habitats using the data 
pooled across three rivers in each season sepa¬ 
rately. The taxa highly associated with habitat 
type were identified using correlation between 
the individual taxon abundance and the canon¬ 
ical axis separating habitats. To calculate the 
mean correlation value for each taxon, these 
correlations were averaged across seasons. To 
examine the community gradients, Principle 
Coordinates Analysis (PCoA) (Gower 1966) 
was employed and indirect ordination plots 
were generated. Both CAP and PCoA were 
conducted using the CAP package developed 
by Anderson (2003). In all cases, the Bray-Cur- 
tis dissimilarity coefficient was used. 

Results 

Physical and chemical parameters associated 
with sampling sites 

The physical and chemical characteristics of 
sites are summarised in Fig. 2. Flow occurred 
in all three rivers during winter and spring 
but there was very little or no flow throughout 
summer and autumn. Bank habitats had higher 
flow than willow root habitats during winter 
and spring when average river flow was high¬ 
est. Percentage CPOM amount in willow root 
habitats were higher than in bare bank habi¬ 
tats throughout all seasons. Percentage FPOM 
showed the same trend within two of the three 
rivers. Slightly lower dissolved oxygen levels 
were recorded in willow root habitats during 
summer and autumn seasons. Slight reduction 
in temperature was detected in willow root hab¬ 
itats during summer. Conductivity and pi I did 


not show significant variation between habitats; 
however, these parameters varied significantly 
across seasons, rivers and sites. 

Effect of habitat on macroinvertebrate taxa 
richness , abundance , biomass and diversity 
There was a significant effect of habitat on spe¬ 
cies richness and log abundance (p>0.05) but 
this was not consistent across the different rivers 
during different seasons of the year (p<0.05 for 
season x river x habitat three-way interaction). 
'Therefore, analysis was carried out separately 
for the different rivers and seasons. Willow root 
mats had higher species richness for most sites 
in all seasons. Total abundance of taxa showed 
a similar trend (Figs 3a and 3b). 

Macroinvertebrate biomass also differed sig¬ 
nificantly (p<0.05) between the two habitats. 
In all seasons, macroinvertebrate biomass was 
higher in willow root habitats (Fig. 3c). 

Shannon Diversity Index (SDI) 

SDI also showed significant differences be¬ 
tween habitats (p<0.05). However, a significant 
river X habitat and season x habitat interaction 
also was observed in the analysis. Therefore, 
analyses were conducted separately for the dif¬ 
ferent rivers during different seasons. Pair-wise 
comparisons indicated an inconsistent pattern 
of SDI between the two habitats during differ¬ 
ent seasons (Fig. 3d). 

Effect of habitat on macroinvertebrate com¬ 
munity assemblages 

Willow root habitats and bare bank macroin¬ 
vertebrate community assemblages showed sig¬ 
nificant differences (p<0.05) in all sites between 
seasons. Correlations of taxa with canonical 
axes showed that eight taxa were associated 
with willow root mats (mean correlation coef¬ 
ficient >0.2). Ten taxa were moderately associ¬ 
ated with root habitat (mean correlation coef¬ 
ficient >0.1) (Table 1). 

In contrast, four taxa were associated with 
bare bank habitats (correlation coefficient >0.2) 
and ten taxa were associated with bare bank 
habitats (correlation coefficient >0.1) (Table 1). 

Discussion 

The differences in physico-chemical parame¬ 
ters associated with willow roots and bare bank 
habitats can be explained by the structural dif¬ 
ferences of the two habitats. Flows recorded 
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conductivity 




DO 



Temperature 





Fig. 2. Mean values of physical and chemical pa¬ 
rameters associated with sampling sites and habitats. 
BC= Birch Creek, JCC= Jim Crow Creek, MR= East 
Moorabool River. 

fl = willow root habitat 

□ = bare bank habitats 

| = stream average flow 
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(a) 


(b) 


Fig. 3. Mean ± SE of (a) species 
richness, (b) log abundance, (c) 
Biomass, and (d) SDl of macroin¬ 
vertebrates, associated with wil¬ 
low root habitats and bare bank 
habitats during four seasons of 
the year. 

^ = willow root habitat 
] = bare bank habitats 


(c) 



Vol 132(4)2015 


101 





























































































































































Research Reports 


s 

JS 

-D 


« 

y 

’£ 

o 

c 

rt 

U 

-C 

(A 

> « 

215 
* « 

2 -c 

o a 

3 U 

° is 
ic « 

So 

3.2 


<-> £ 
3 3 

213 


d 

B 


cu t: 
< o 

u g 

§ * 
2 o 
<0 3S 

J2> 
I O 

-D c 

Q> O 
S? rt 

-3 CU 

*73 O 

g « 

cn 2 
« « 
.2 2 
15 ^ 

J2 G 

(A 


$% 
<- o 

•S CX 


e 

§3 
« S 

■g£ 

4/5 r5 

VO O 

n 3 
« c 
* O 

£ >3 


DC 

.5 

Dm 

C /5 


Jcr^t^-\omm-r-rror<~i-^ on o\ oo ao oo 

— OOOOO 

dddddddddddddddddddodddddddd 


On ’H On (N M V) rO 

-T-rrM-rTrpr^fN 

dddddddd 


ON 00 
O cm 

o d 


in 03 - OO (N cc rn 

O _ _M O —' ~ ~ 

d> d> <£ d> d> o o 


r-. o\ 

CM CM 

d d 


o o o o o 


>—< —■ m in \c m cm cm ^ 

—' CM CM —< —< —' —« —< CM —< CM 

ddddddddoddo 


-H Ln \C -H On 't rr n t 

rc>cMrOrc>F—(CMCMCM- hCMCM 

dddddddddod 


OO rO MD oc 
rO O CM O 

dodo 


o o — 


oooooooo 


Cm k. 
v> p 

d, 2 d 

x S ^ £? S 

5b § S d « 

<u "S "rs 

•§ £? f a 
£|;ts | 
o' 3 -5 <£ 6 


CM MO ON 
CM O O 

dddddddddd 


.3 

IS 


3 


O 

o 


d, -s: 

vT Om 
Q 3 

,= 5 t 


cu 


3 

U 

<3 

S; 




g>3 

B H 


a> O 

U U 


4 » 

Q> 

«9 

0 / 

« 

T 3 

<u 

« 

a> 

C 3 

a» 

T 3 

-o 

•5 

3 

3 

C 3 

3 

IS 


Tj 

3 

U 

CJ 

*c 

o 

c 

£ 

p 

'A 

-a 

c 

cl 

£ 

n 

« 

i 

rt 

i— 

« 

U 5 

Dm 

s 


CU 

X 

CU 

ua 


SS CU <U 

-2 « « 

2 2 u 3 

.a 3c -o 32 
£0.^0, 

§• a .§ S 

S. E £ S 

3 iTu 2T 


» a 

-S 3 
U 

r 1 * tn 


<U 

C3 

3 

2 

cu 


Q> 

CJ 

*3 


C5 

CU "O 
« — 


-a 

^ 2 *r £ 

V 5 < D- 


wrj V 

Cu J2 •- « 'a 

O T3 GS T3 C 

‘ • «J 3 ’O O 

^ S .b ~ 


11 
.3 cu U> 3? 
73 t/5 U1 ^ 


CU 

« 

’k 

cu 

u 

o 


CU 
ai 


rt T3 

3 c 
-o o 
O ‘C 

2P 
o a 

£ "O 
3 O 
« cu Q. 

>. So 

• O & 

Dm ^ 


44 a< 

^ CJ 

Dm O 
W *4 

E O 
£ 

^ X 


102 


The Victorian Naturalist 








Elmidae (larvae) 0 03 0.23 0.04 0.01 0.08 

Calamoceratidae 0.1 0.05 - 0.06 0.07 
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Leptophlebidae sp. Wundacaenis sp. -0.06 0.1 - - 0.02 

Glossiphoniidae 0.02 0.1 0 -0.03 0.02 





Bare banks Sphaeriidae Sphaeriumsp. - 0.03 - 0.44 - 0.23 - 0.29 - 0.25 

Psephenidae Sclerocyphott $p. - 0.37 - 0.03 - - 0.27 - 0.22 

Empididae - 0.22 - - - - 0.22 
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Sailidae 0.13 -0.18 - - -0.02 

Lymnaeidae Glyptophysa smooth -0.04 0.05 -0.06 -0.04 -0.02 

Lestidae Austrolestes a nidis -0.02 - - - -0.02 
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in willow root habitats during winter and spring 
were comparatively lower than in bare bank habi¬ 
tats in most sites sampled, suggesting that willow 
roots would act as mechanical barriers and filter 
suspended material. Higher amounts of CPOM 
and FPOM occurred in willow root habitats com¬ 
pared to bare bank habitats, even in winter when 
willows had lost their leaves, further confirming 
that willow root mats filter suspended material. 
The relatively reduced temperature levels recorded 
in willow habitats in summer would be due to the 
shading effects of the canopy (Pidgeon 1978; Glova 
and Sagar 1994). r Ihe lower dissolved oxygen lev¬ 
els of willow root habitats during summer and 
autumn can be related to the increased levels of 
organic matter coupled with reduced stream flow 
in these seasons. Increased community respiration 
coupled with reduced dissolved oxygen levels as¬ 
sociated with stagnant water created an oxygen 
gradient between the two habitats. This also could 
have been affected by the heavy shading of willow 
habitats during summer and autumn as willows 
would have had a full canopy that would have re¬ 
duced the photosynthetic algal growth under wil¬ 
lows compared to open bank habitats (Pidgeon 
1978; Glova and Sagar 1994). 

The interaction of season, river and habitat on spe¬ 
cies richness and abundance indicates that the effect 
of habitat is not consistent in different rivers during 
the different seasons. These differences may be due 
to inherent heterogeneity of available resources or 
patchiness of each river, and governed by their flow 
regimes. Willow root habitats are more favourable 
for macroinvertebrates in terms of food availability 
and refuges compared to mineral substrate. Simi¬ 
larly wallow root mats act as mechanical barriers, 
which prevent the effects of high flow and stabilise 
structures for their attachment compared to the 
more unstable mineral substrates. Total abundance 
of taxa showed inconsistent patterns from site to 
site, probably due to the different compositions 
of taxa in different habitats. Some micro levels of 
physico-chemical parameters favour the abundance 
of particular taxa, which leads to their dominance, 
creating more variability in total abundance between 
sites, Besley (1992) compared riparian root habitats 
of River Red Gum Eucalyptus cam aid i den sis. River 
Oak Casuarina cunninghamiana* and White Willow 
Salix alba with bare bank habitats and showed that 
riparian root habitats supported significantly more 
species than bare bank habitats during autumn and 
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winter. Results of the present study were consist¬ 
ent with these findings. 

Total biomass of taxa was consistent in willow 
root habitats and bare bank habitats during all 
seasons and rivers, although higher biomass of 
taxa was observed in willow root habitats in all 
seasons. Glova and Sagar (1994) showed higher 
species richness, total abundance and biomass 
of benthic invertebrates associated with willow 
than non-willow sections, consistent with our 
results. The increase of fish abundance associ¬ 
ated with willow habitats in the previous study 
may be attributed to the higher food levels as¬ 
sociated with willow habitats 

SD1 is more affected by the presence and abun¬ 
dance of taxa. Lowered levels* of SD1 in willow 
root habitats in some sites indicate that few taxa 
dominate that habitat. This suggests that wil¬ 
low root mats favour particular taxa and their 
abundance. Organic matter enrichment in wil¬ 
low root habitats favours pollution tolerant taxa 
by replacing more sensitive taxa groups (Suter 
1990). Fauna of willow root mats was dominat¬ 
ed by deposit feeding collectors such as Oligo- 
chaetes and Amphipoda (Latta 1974). Plecop- 
tera, Trichoptera and Ephemeroptera, which 
respire with gills or direct cuticular exchange, 
can be particularly susceptible to lowered dis¬ 
solved oxygen levels (Dallas and Day 1993) and 
the taxa which were missing from willow sites 
in summer were consistent with this pattern. 

Most taxa showed higher association with wil¬ 
low root habitats during autumn when most of 
the allochthonous organic matter in the form of 
food from willow reaches the stream. Ihe dif¬ 
ferences in community assemblages observed 
in the two habitats can be related to many fac¬ 
tors. Amphipoda (Paramelitidae), a relatively 
large group of shredders, was the dominant 
group in willow root habitats in all seasons and 
were highly associated with willow root habitats 
during autumn, suggesting an association with 
CPOM content. Ceinidae were highly associat¬ 
ed with willow root habitats during most of the 
seasons suggesting that they are dependent on 
FPOM or decaying root mats or have a structur¬ 
al preference for this habitat. The higher associ¬ 
ation of more organic pollution tolerating taxa, 
such as Potamopyrgus sp., Cura sp., Megadrilii, 
Physa acuta and shredders such as Amphipoda, 
with willow root habitats during autumn in this 


study suggests that they prefer organic matter 
enrichment available from willows during au¬ 
tumn. The higher association of ambush preda¬ 
tors, such as Coenagrionidae, with willow root 
mats suggests that predators also benefit from 
the complex structures of root mats. 

Simuliuni sp. and Leptoperla sp. showed mod¬ 
erate mean correlations during winter and 
spring. Sinntttum sp. uses cephalic fans for cap¬ 
turing food items and may prefer the higher 
flows of winter, which facilitate its food captur¬ 
ing habit. Leptoperla sp. was also highly associ¬ 
ated with willow' roots during w inter and spring. 
Its preference for this habitat may be related to 
well-oxygenated fine particulate organic mat¬ 
ter in the rich willow root environment. Some 
grazers, such as Gyrauls sp., Glyptophysa sp., 
Potamophyrgus sp., were associated with willow 
root habitats during winter and spring, benefit¬ 
ing from periphyton and biofilm associated with 
willow root mats under the open canopy of wil¬ 
low's and from the stable substrate of willow root 
mats, which provide refuge during high flow. 

In contrast, Sphcarium sp., Psephanidae, Tu- 
bificidae, Corixidae, Phreodriliidae, Ancylidae, 
Pisidiutn sp., Tipuliidae, Capilariventridae, 
Dytiscidae Lumbriculidae, and Leptophleb- 
biidae w'ere highly associated with open bank 
habitats. Their association with mineral sub¬ 
strate may be related to their food acquiring 
behaviour. Most of them are filter collectors, 
which benefit through such mineral substrate 
and associated suspended organic matter. Some 
predators such as Corrixidae and Dytiscidae 
are also associated with bare bank habitats. 

Ihis evidence suggests that willow' root mats 
are productive habitats compared to open bank 
mineral substrates. This habitat supported 
distinct macroinvertebrate communities and 
higher species richness and biomass compared 
to bare bank habitats. This suggests that remov¬ 
al of willow's from streams where other ripar¬ 
ian cover is lacking can have detrimental effects 
on stream macroinvertebrate communities. 
Boulton and Lloyd (1991) suggest the need of 
alternative mechanisms to improve hahitat het¬ 
erogeneity where willow removal takes place. 
Further research is needed to understand alter¬ 
native methods, such as macrophytes, snags or 
fast growing vegetation cover, e.g. reeds, would 
all be effective in improving habitat heteroge- 
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neity and thus biodiversity in streams where 
willow removal is taking place. 
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